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Abstract 

We develop a variational formalism in order to study the structure 
of low energy spectra of frustrated quantum spin systems. It is first 
applied to trial wavefunctions of ladders with one spin 1/2 on each site. 
We determine energy minima of these states. The variational ground 
state shows a finite energy gap with respect to the energies of states 
which span the Hilbert space and are orthogonal to it. This is the case 
for any size of the system. Under some justifiable approximations the 
argument can be extended to even-legged ladder systems in 2d and 
higher dimensional spaces. The Hamiltonian can contain spin-spin 
coupling interactions of any range. For specific values of the coupling 
strengths level degeneracies can occur. 
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1 Introduction. 



The low energy spectra of magnetic materials can be described by quantum 
spin models characterized by their geometric structure, the form of the in- 
teractions between the spins located on sites, the strength and the range of 
the couplings which fix the intensity of the interactions. The knowledge of 
the spectral properties of these systems is considered to be of prime impor- 
tance for the understanding of specific phenomena such as superconductivity 
at high T c . Investigations on this subject which go on for many years have 
concentrated in the recent past on specific structures, in particular the pres- 
ence of ladders and stripes in 2d superconducting material, see f .i. [TJ |21 El H] ■ 

It has been observed that systems which expectedly behave like 2-leg lad- 
ders show a gap between the ground state and the first excited state [3] 
whereas those which behave like 3-leg ladders possess a continuous spec- 
trum jB]. Theoretical investigations have been developed in order to study 
the properties of the ground state and low energy states with a particular in- 
terest for the existence or absence of such a gap. In Id systems Lieb, Schultz 
and Mattis (LSM) [7] established that spin-1/2 systems which are invariant 
under translations and rotations either show no finite energy gap between 
the ground state and the first excited state, or the ground state is degenerate 
and spontaneouly breaks the translational invariance creating a dimerized 
structure. In the first case the energy interval between the ground state and 
the first excited state AE goes to zero as 1/L when the length L of the chain 
tends to oo. A large amount of work on spectral properties of quantum spin 
systems followed and can be found in the litterature from the late 60s to the 
late 80s, see f. i. |H1 El UQI EED EES • 

An important step was performed by Haldane who conjectured that the 
spectra of Heisenberg antiferromagnetic chains are gapless if the spins are 
half-integers and show a gap if the spins are integers ^1E]- Later Shelton 
et al. |T5] found that weakly coupled isotropic spin-1/2 antiferromagnetic 
Heisenberg chains always show a spectral gap. In an application of the non- 
linear sigma-model Sierra jT^j tried to extend the Haldane conjecture to lad- 
ders. He showed that in the framework of this model spin-1/2 even ladders 
exhibit a gapped spectrum whereas odd ladders show a gapless spectrum. 
Further tests of Haldane's conjecture have been made, see f.i. jT71 [THj and 
refs. therein. Oshikawa extended the LSM theorem to quantum spin systems 
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with a conserved particle number on a periodic lattice in arbitrary dimen- 
sions and showed that an excitation gap opens when the number of particles 
per unit cell of the ground state is an integer [T^]. The theorem was also 
found to be verified in chiral liquid spin systems in higher dimensions [2*U] . 
There it was however observed that the conclusions drawn in ref. [01] may 
not be universally verified. More recently Hastings |2I] showed that AE 
goes to zero like InL/L for spin- 1/2 systems in higher dimensions when the 
Hamiltonian of the system is SU(2) invariant and sites are occupied by an 
odd number of spin- 1/2 particles. 

Further information has been obtained by means of specific models, in 
particular quantum spin systems, ladders and 2d networks, see in particular 
[22 EII1 EH EI3 Oil • These studies deal with specific systems using different 
formal tools. They show non trivial low energy spectral properties and very 
rich phase diagrams. However the low energy properties of 2d and higher 
quantum spin systems remain an open question in a large number of cases. 

In the present work we propose a further investigation concerning the 
existence of a finite energy gap AE between the ground state and the low 
energy states of frustrated quantum spin systems. More specifically we con- 
sider spin 1/2 ladders and extend the analysis to even- legged 2d systems, as 
well as systems with an even number of chains in an arbitrary number of 
space dimensions. We show that the results are valid for any range of the 
spin-spin interaction. 

The investigation relies on a variational approach. We consider first the 
case of spin-1/2 ladders. The central concept which governs the behaviour 
of the system is SO(4)-invariance of the Hamiltonian. We introduce a trial 
wavefunction in order to determine the energy extremum of the system and 
analyze the nature of this extremum. If this state corresponds to a mini- 
mum in energy we show that the states orthogonal to it are characterized 
by an energy gap with respect to this state except for specific values of the 
coupling constants which enter the Hamiltonian leading to degeneracies. We 
then show how the results can be extended to 2d and higher dimensional 
spaces for systems with an even number of legs. Finally we comment and 
discuss the outcome of the present investigations and draw conclusions. 
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2 The model. 



We consider spin-1/2 ladders (23121] described by Hamiltonians of the fol- 
lowing type 

N 

H (s,s) _ J t y^ ShSh + J l s^S^ + Ji } j S i2 S i2 + Ji c 2j Silvia (1) 
i=l <«i> <«j> <*?> 

+ ^2c 2j S *2 S .h 
<ij> 

where the indices 1 or 2 label the spin 1/2 vector operators Si k acting on the 
sites % on both ends of a rung, in the second and third term i and j label 
nearest neighbours along the legs of the ladder. The fourth and fifth term 
correspond to diagonal interactions between nearest sites located on differ- 
ent legs. N is the number of sites on a ladder and the coupling strengths 
J t , Ji, Jic, Jic are positive. In the sequel we fix J\ c = J 2c = J c . The ladder is 
represented in Fig. 1. 

The Hamiltonian (1) possesses SO (4) symmetry by construction. By 
means of a spin rotation [3S] 

Si^^Si + Ri) . (2) 

Si 2 = 7}(Si ~ Ri) • (3) 
it can be expressed in the form 

T - 

H (S,R) = ± _ R 2 } + Ji ^ Sfi + J 2 J2 (4) 

i=l <ij> <ij> 

where the components S^ + \ \ and R\ , R\ \R^ of the vector op- 
erators Si and Ri are the SO (A) group generators which can be written as 



s\ +) = V2pr, (11)(10) + x i (10)(1 " 1) ) = s\-> 



s [z) = x (11)(11) - x (1_1){1_1) 
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Figure 1: Top: the original spin ladder. The coupling strengths are indicated 
as given in the text. Bottom: The ladder in the SO(4) representation. See 
the text. 

R { : Z) = _(X (10)( ° 0) + X (0 ° )(10) ) 

where 

xf M ° ){s ' Ms,) = \SM s ) l . t (S'M s ,\ 
and the states \SM S ) are generated by the coupling of spin 1/2 states 
\SM s ) i= J] (1/2 mi 1/2 m 2 |10)|l/2 7^)^1/2 m 2 ), 

mi,m2 

to two particle singlet (S = 0, M s = 0) and triplet (S = 1,M S = 0, +1, -1) 
states. 

In the new representation the coupling stengths are J\ = (Ji + J c )/2 and 
■h = (Ji — Jc)/2- Ji is always positive, the sign of J 2 depends on the relative 
strength of the couplings between spins on sites located on the same leg (J;) 
and different legs (J c ). The 5*0(4) symmetry is shown to be of particular 
interest for the characterization of the spectral properties of even-legged spin 
systems. This point appears clearly below. 
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3 Variational ansatz. 



3.1 N-body wavefunction. 

Variational techniques have been applied for a long time in order to investi- 
gate the properties of quantum spin systems, see f. i. refs. [TTJ I3UJ EH E2J 
HUH IHH HHH HUE EH] ■ In analogy with the phenomenological BCS approach and 
the approach of refs. (3^1 EH |2H I2S] we introduce a trial wavefunction which 
takes a product form of two spin- 1/2 states 

TV 

i«o w >=ni*i ro >- < 5 > 

i=i 

Here l^ "*) is the wavefunction at the location i which is chosen as a linear 
combination of singlet and triplet states \SM s )i 

|^ (0) ) =a i |00) i + /3 ff |l-l) i + /3 io |10) i + Ai|l + l)i • (6) 
with the normalization 

l«d 2 + IAtT + lAol 2 + lAil 2 = i - (7) 

The time-reversal symmetry of the Hamiltonian allows to introduce |/%| 2 = 
|/3ji| 2 . A further constraint is imposed such that fin = /% = K(3 i0 where k is 
fixed as the ratio 

(1/2 1/2 1/2 1/2|1 1) /- 
(1/2 1/2 1/2 - 1/2|1 0) ' 1 J 

where the Clebsh-Gordan coefficients weigh the coupling of the two electron 
state \l/2m\\ l/2rri2)i to the states |ll)j and 1 10) ^ respectively. Each site is 
occupied by a single spin. The constraint takes care of the relative weight of 
the |1 1) and the |1 0) two-spin states. Defining 7, = (1 + 2K 2 ) 1 / 2 /3 i0 leads to 
the normalization relation 

\ ai \ 2 + | 7i | 2 = 1 . (9) 
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The imposed constraint fixes (3n = k^/ (1 + 2k 2 ) 1 / 2 = /%. The {a^} and 
{7;} are the variational quantities. 

Using the expressions of the operators Si and Ri defined above the en- 
ergy = \H^ S ' R ^\^^) associate to the trial wavefunction |$q ) * s 
obtained as 

N N 

E° N = 2E T + E s ^( 3 l«i! 2 + M 2 ) + Einti + E int2 . (10) 

i=l i=l 

with E T = J t /4, E s = -J t /4, 



and 



;i + 2k 2 ) 2 

v ' <ij> 



Eint2 = (1 + 2 K 2) E^WW ^ 

+ 7iril 2 y J2( Im ( a ^ Im ( a ^ 

^ ' <ij> 



3.2 Extrema of the energy 

In order to fix the energy spectrum of the system the complex amplitudes a.{ 
and ji are varied under the normalization constraint Eq. (JOJ) 

' + = • (13) 



da* da* c?7i 

We introduce the parametrization 

^ = exp(i<f> u )sin(9i/2) . (14) 

and 

7i = exp(i4>2i)cos(8i/2) . (15) 
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Applying Eq. (JT3J) to E° N given by Eq. (jlUJI and using Eq. (|14|) and Eq. (|15j) 

leads to 

eiexp(i((j)u - <fai)8in(9i) + (cos 2 (^/2) - sin 2 (0 i /2))A i = . (16) 

where 

e^^-^-^J^L,^!^. (17) 
The energy q is a negative quantity, A, is complex and given by 

1 + 2k 2 



A * = - 2 « 2 K*7i + (1 + 2«> i7 *) • (18) 



The symbol < j > indicates the nearest neighbour of i, in practice j — i + 1. 
Aj can be written as 

Ai = lAilexp^^M - 2i )) . (19) 
The solution of Eq. (|T6*j) takes the form 

eiSinOi + | Ai| cosdi = (20) 

with 

sm0< = -|A. 4 |/^ (21) 

and 

cos9 i = e i /E i . (22) 

Here £ 2 = e 2 + |Ai| 2 and 

l^| 2 = ^(l + |)- (23) 

N| 2 = ^(l-|)- (24) 

A, = A[ + <A* . (25) 

A r = ^ T^-. (26) 

V A| 

~~ (1 + 2^)^^- (27) 



The energy Ei can in principle be positive or negative. The expression of 
6j in Eq. (JT7j) shows that the variational wavefunction |$g ) postulated in 
Eq. (0) entangles nearest neighbouring site state wavefunctions. 

The Hamiltonian of the spin ladder is fixed by the three coupling con- 
stants J t , J\ and Ji. A close inspection of the expressions of (Eq. (fT7jl ). 
Aj (Eq. (|T%j)) shows that the amplitude a« and hence ji of each site can be 
determined self-consistently if the system is closed in such a way that site 
iV + 1 is identified with site 1. 



3.3 Nature of the extrema. 

The nature of the extremum is given by the sign of the second derivative 
d 2 E^ of the energy E^ with respect to under the constraint fixed by 
Eq. Q and the definition of 7«. The explicit calculation leads to 

4 |7i| 2 It* I i 1 + 2k ) ^ 

The sign of d 2 E^ depends on the relative strengths of J t , J\ which are positive 
quantities and the ratio of the amplitudes and |7j| 2 . In fact 

ff = & (29) 

which is larger or smaller than 1 depending on the sign of E^. 

d 2 E% is positive if Jt is small enough so that the first term in Eq. 
gets smaller than the second term which is strictly positive. 

More precisely, for fixed J t and J\ d 2 E° N is positive if |aj| 2 /|7i| 2 gets larger 
than a minimum value, 

|a - |2/N2 a 1 + (30) 
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where 



This may be realized if J\ is large and (or) |Aj| is small enough which may 
be the case for small |J 2 | = \{J\ — J c )\/%- 



4 Energy of the variational state. Nature of 
the spectrum. 

4.1 Unprojected wavefunction. 

Using the explicit expressions of the amplitudes {ctj} and {7j} leads to the 
expression of the variational energy 

E% = (E T + 2E s )N + (E T -Es)J2^r (31) 

v ' i <j> J 

(1 + 2K 2) 2^e.Z^ E . (i + 2 «2) ^e.Z^ Ej 

We want to analyze the properties of the energy spectrum, in particular 
to see whether it is discrete or continuous. 

Since the Hamiltonian of the system is real symmetric the space of states 
is spanned by a complete set of states which are orthogonal to each other. 



N 

|*W) = JJ|„<* 



i=l 



Orthogonality is realized if \&^} differs from |$o^) by a ^ least one state 
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|*f> = ai fc) |00> i + 7 ? } |10>,. (32) 

orthogonal to 

|vl/f ) ) = ai (0) |00) l + 7, (0) |10),. (33) 
where af 1 and 7^ stands for the former notation ai and 7$ of section 3. 

If a\ = —Ji and = at$ * the wavefunction in which is 
replaced by ) is then orthogonal to it by construction. Using the 
energy ^u-i) corresponding to the location of site (i — 1) of the state N ') 

orthogonal to the corresponding state l^f^) in |$q 7V ' ) ) is then given by E<j> 
corresponds in fact to the site i, see below Eq. fITHj) 



^, = (^-^)-(i^i#r (34) 

where |7,- | 2 = (a- '*! 2 . Eq. (34) shows that in general 

4l> * «&) • PS) 
Using Eq. (JTBJ) it is easy to see that 

AS = -AS • (36) 

and hence 

4-D * <d • ( 37 ) 



Consequently the states l^^) and l^f^) are non degenerate states since 

e (fc) £ (o) 

= {Et-EsX-™-™-) (38) 



2 w (°) 

j \r A r ( k ) A r (°) 



(1 + 2^^ j^' 

2« 2 J 2 Aj Ag Ag_ 
(l + 2/t 2 )£ i V fc) J 

v ' 1 (i-1) 
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is different from whatever the number of sites except if (A^j-J = |A| t -Ji| = 
which may happen if J 2 = i. e. Ji = J c . 

The physical states which are orthogonal to the variational state are gen- 
erally linear combinations of the orthogonal states Since these states 
are non-degenerate with the variational ground state this property will be 
shared by the physical states in the general case, i.e. there is a finite gap 
between the variational state and the other physical states. As already men- 
tioned above degeneracy in energy may occur in special cases, f.i. if J 2 = 0. 



4.2 Projected wavef unction. 

Since the Hamiltonian H^ S ' R ^> commutes with the total spin projection M tot 
the actual wavefunction has to be projected on a state of fixed total spin 
projection 

1 r ir N A 

= ±- / rf^-S.^nE^I^^ . (39) 

2irJ -* Lit! 

where r\\ = a u r\\ = /%, rj l 3 = (3 i0 , r\\ = (3 n . 

The energy corresponding to the projected state V\<&^) is then obtained 
from the hermitic matrix element 

E% = h{$W\H^V\<t>l N) ) + (& N) \VH^\<t>l N) )} . (40) 
In the case where M tot = one gets 



(*W|fl<wip|.W) = J2-^^TT £ 4,.,A(h.-,Q ■ (4i) 

where c\ , is the combinatorial coefficient associated with the decom- 
position of ^(Sf )'«, I = 2p = ^ h and 
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TV 4 

1=1 fc=l 

where {M^i} are the projections of the states k on the sites %. 

Using the explicit expressions of the {ttj}, {7^} and {A«} (eqs. (23) to 
(27)) £%{li, .., l n ) can be put in the form 

■■, U = + - 2 (1 + (-l)'l H\ 2 (43) 



A 2 

i 

"A 4 



<ij> 

+ J E ^ i) a *i a ]lil3 - VijQUiUjltl* + B(j, 1)^0*^1 + j)«*«i7i7j-] 

<ij> 



where A = (1 + 2/t 2 ) 1 / 2 and 

= « 2 [l + (-l) ,i+ ^ + (-l)' i+1 + (-l)^ +1 ] (44) 
+26 li0 6 lj0 + S h0 (l + (-l)' i +2(-l) i< ) 

fl(<j') = M 1 + (- 1 )''+*i.o] • ( 45 ) 

and 

C(i,j) = + + (46) 



One can then construct projected states {Pl^j^)} which are orthogonal 
to V\Q$) in the same way as in section 4a) above and find out that the 
energy difference AS = £^ — £q and the corresponding AEp = E P — E P 
are different from zero except maybe if J 2 — 0. 



13 



5 Generalization to even- legged spin 1/2 sys- 
tems in 2d with any type (short or long) 
range interaction. 

It is straightforward to see that the present considerations can be applied 
to spin ladders which do not only contain interactions between neighbouring 
sites, but also interactions at any distance from each other. Such contribu- 
tions introduce a quantitative but no qualitative change in the gap equations 
Eqs. (J2S1), (EED and Q27j). 

More generally, the results can be extended to 2d frustrated systems with 
an even number of ladder legs showing the same type of couplings between 
sites as those introduced on the ladder. They are obtained by adding ex- 
plicit couplings between sites located on the neighbouring legs which belong 
to neighbouring ladders. If L = M * N where M is the number of parallel 
ladders and iV the length of the ladders, the Hamiltonian in the spin rotated 
representation (Eqs. (2) and (3)) can be put in the form 

H (s,R) = Jj_ ST {S * - Rl) + S l S J + J 2 J2 R i R i + H[mt) ■ (47) 

i=l <ij> <ij> 

with = if M 1 ) + and 

jr L-N 

jj(tntl) _ _j^t (SjSj + N + SiRi + N — RiSi + N — RiRi+N) ■ (48) 
i=l 



J: 



jj(int.2) int 



'L-N-l 



l+N 



i+N+1 



(49) 



i=i 



—RiSi + N+l — RiRi+l+N) 
L-N 

+ {£>i-lSi+N + Si-lRi+N ~ Ri-lSi+N ~ Ri-lRi+N) 



i=2 



J\ nt is a coupling strength between sites in neighbouring chains belonging to 
neighbouring ladders along rungs and Jf t acts between the same elements 
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along a diagonal linking neighbouring site (see Fig. 2). 

The Hamiltonian Eq. (JTTj) to Eq. describes a unique chain which is 
constructed in such a way that the end of each chain m (m = 1, ...M) corre- 
sponding to the Hamiltonian H^ S ' RS) Eq. (JJj) is connected to the beginning of 
the next one, see Fig. (3). 

This introduces spurious terms S kN S( k+ x) N , RkNR(k+i)N, RkNS( k+1 ) N , 
SkNR{k+i)N, k = 1, M — 1 which couple the border sites of the consecutive 
spin chains. The contributions of their matrix elements should be subtracted 
from the expression of the energy after the variation. The number of terms 
of this type is however small compared to the total number of states if M 
and N are large ("surface" contributions). Neglecting them in the limit 
L = M * N — > oo one can go through the calculations as before, obtain the 
same type of equations and draw the same type of conclusions as before. 
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i = (M - 1) N+ 1 



Figure 3: The 2d system of fig. 2 in the SO (4) representation. The long 
dashed links are the spurious contributions. See explanations and discussion 
in the text. 



More explicitly the expression of the energies is given by 

v ' <j> 

Jl jr _ 9^.2 Td 

and the complex gap A, which is now the sum of three terms 

\ = A u + A 2i + A 3i . (51) 

■h 

<j> 

^ = ~ 4(1 ^ 2) - 2 « 2 )«i7i + (! + 2 "V/>?) • (53) 

(i) 

A 3* = ~ 2(1 f^ 2) E((! - 2 *>^ + (1 + 2^7*) • (54) 

< j > (= i + 1) stands for the site which is the right-nearest neighbour 
along a chain , (J) corresponds to the nearest-neighbour of a site % along a 



where 



VH = ~ (1+ J 2«») " 2K2) ^ 7j + (1 + ■ (52) 
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rung between chains belonging to nearest-neighbour ladders and [j] to the 
cross-diagonal nearest-neighbour of site % located on the chain of the nearest- 
neighbour ladder. 



The expression of the variational energy is formally the same as the one 
of section 4 except for the additional sum of terms corresponding to the 
couplings JJ mt and Jf nt . It is given by 

U) bl 

The nature of the extremum (maximum or minimum) depends on J t , J\ and 
the additional coupling strengths J\ nt and Jf nt . d 2 E N is positive if |J 2 | = 
\Ji — ./cj/2, J r int and Jf nt are such that 

N7N 2 > , , _ 3 

Jt ^r\r*, ui, u int , u int 

where 

2k 



1 + -j-G(k, Ji, J\ nV Jf r 

G(K, Ji, Jl nt , Jfj = (1 ^ 2)2 [4J 1 M + Jlnt E l^'l 2 + 2J "nt E 



(55) 



Hence d 2 E N is positive if Ji, Jf nt , J^ t get large and J t and (or) the {Aj}s 
remain relatively small. 

The discussion about the structure of the spectrum can be taken from 
section 4. In general this spectrum will show a gap between the ground state 
and excited states, except for specific values of specific coupling strengths. In 
fact degeneracy of states appears when the gaps A ni (n = 1, 2, 3) disappear, 
which happens when J 2 , J r int and Jf nt are equal to zero. 

The present results can evidently be extended to the case where one works 
with projected wavefunctions along the same lines as it has been shown in 
section 4. 
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6 Remarks, discussion. 



The present derivations are subject to comments: 

• The spin Hamiltonians H which were introduced commute with also 
with the total spin S to t- Hence the corresponding quantum numbers 
are conserved and characterize the states of the physical spectrum 
l^o^ {Stot) M tot )). The trial wavefunction |$q ) should be projected 
on a state with fixed S tot . This changes quantitatively the energy E° N 
and Ep. It should however not qualitatively affect the calculated energy 
differences and the existence of a gap between the variational ground 
state and other states of the infinite system. 



• In practice (3n = fin may be treated as independent variational param- 
eters if the constraint introduced in section 3a) Eq. (|SJ) is not taken into 
account. The corresponding component may be kept in the definition 
of |^ 0) ), Eq. ®. 



• It is easy to see that the present formalism can be extended to any 
even-legged system of spins in any space dimension D = nd with the 
proviso that, as in the 2d case, the contribution of terms coming from 
the border sites of the finite systems and which should not be present 
can be neglected when applying the variational procedure. This is 
realized if the system gets large. Since the algebraic expressions are 
rather cumbersome we do not write them out explicitly. 

• The "linking procedure" used in Section 5 in order to generate a con- 
tinuous linear system leads a posteriori to a symmetry constraint of 2d 
and systems of higher dimensionality which by construction possess the 
5*0(4) symmetry introduced in section 2. This is a constraint which 
may limit the range of the variational procedure in the sense that other 
symmetries could take over at different space dimensionalities. 
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7 Degeneracies. 



7.1 Behaviour of e$ and Aj at J2 = 0. 

As observed above in sections 4 and 5 the variational state gets degenerate 
with states orthogonal to it when the gaps {Aj} vanish. This is the case 
when J 2 = (section 4), J 2 = Jf nt = J^ t = (section 5). It is of interest 
to analyze more precisely the property of this point. In the sequel we con- 
centrate on the case of a ladder system. The generalization to 2d and higher 
dimensions is straightforward. 

It is easy to determine the behaviour of Aj and e$ in the neighbourhood 
of J 2 = 0. The derivative of each gap with respect to J 2 reads 



dA t 1 A m (J 2 ) J 2 . A[(J 2 



(56) 



dJ 2 1 + 2k 2 E i+1 (J 2 ) 1 + 2k 2 y E i+1 (J 2 ) 

a , T , e i+1 (J 2 )e' i+1 {J 2 ) + A m (J 2 )A m (J 2 ) A 

" Am(j2) ^5 } 

where A;( J 2 ) = dA l {J 2 )/dJ 2 and J 2 ) = de^) / d.J 2 . A\(J 2 ) = for J 2 = 
and is a continuous function of J 2 in the neighbourhood of this point. 

Similarly e^(J 2 ) is a continuous function of J 2 in the neighbourhood of 
this point. 



7.2 Eigenstates at J2 = 0. 

We show that the states {1$^)}, k — and k ^ are eigenstates of the 
Hamiltonian H^ S ' R ^ defined in Eq.(4) when J 2 = 0. 

Indeed the part of the Hamiltonian proportional to J t is diagonal in the 
basis of states {|<E>j^)}, the third term does not contribute and the second 
term is a sum of terms of the form 

(^\S,S,\^) = * 7 f' 7 f> 7 f > . (57) 

Using the expression of 7^ in Eq. (23) and the negative root solution of 
Ei, Et = -(e 2 + I Ail 2 ) 1 / 2 leads to (^\SiSj\^) = for any k and k', i 
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and j. Hence the states {|3>jj. }} are eigenstates of H^ S > R ^ orthogonal to each 
other. 



7.3 Behaviour of the energy interval between the states 
in the neighbourhood of J2 = 0. 

The energy intervals between the physical states which are linear combina- 
tions of states {|^ )} vary continuously with J 2 in the neighbourhood of 
J 2 = 0. 

Using the expression given in Eq. (jHEJ) for a state |3>jj. ) which differs 
from the state |3>j^) by n consecutive site wavefunctions \^\ k ^), (f.i. I = 
N — n + 1, N) the energy differences between these states are given by 



A£f' } = EP-E^ = (E T -E s )J2(el k) -el k,) ) (58) 

1 

2K 2 Jl V"^r/ 1 w (fc) (fc')\ 



;1 + 2k 2 ) 2 t-f 



. / (fe) . (fc) / (fc') . (Jfc')s i (fc) (fc) (A') (fc')l 

I J2 \xr (xr{k) _ Mk') x Mk) ( Mk) _ Mk')^ 

TT^2kF) N ~ n+1 N ~ n + 1 ' s - 1 ^ >i 



2k 1 J, 



i N~n\°N-n+l °N-n+l)^^°l V°l+1 



(1 + 2K 2 ) 

where ef ] = e {k) / E[ k \ 6^ = A\ {k) /E[ k) and 8\ {k) = ^ {k) / E \ k) depend im- 
plicitly on J 2 . 

AE^ k ^ goes continuously to zero for any {k, k'} when J 2 goes to zero 
which shows that the degeneracy of states for this particular value of the 
coupling J 2 . The same is true for the physical states which go continuously 
over to the states )} when J 2 goes to zero. This is also the case when 

the states \^f\ ) which differ in different states are non consecutive as 

considered above. 
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Close to J 2 = E$ and hence AE^ k ^ behave like cast * Jf to lowest 
order in J 2 and E%\j 2 ) ~ E<£\-J 2 ) 

Finally the same reasoning is valid in the projected framework dealing 
with Pl^) andP|$jv°). 

8 Generalization to a correlated variational 
wave function. 

8.1 Correlated wavefunction. 

The variational wavefunction which was postulated in Eq. (j5J) has a mean- 
field structure which includes implicit correlations between site states through 
Eq. (J2*U|) to Eq. (J27j) but no explicit correlations. In order to implement such 
correlations we go over to a pseudo-fermion representation of the operators 

s\ +) = b+d, + dfa 

S\- ] = c+di + dfk 
S\ z) = bfbt + c+ Ci + dfd t 
R\ +) = bjdi + aid 
Ri = a\bi + cfdi 
Ri = cfcii + afdi 

where af,bf,cf,df are anticommuting fermion operators which generate 
singlet and triplet states 
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a+|0) i= \S = 0,M = 0)i 

bf\0) l = \S = l,M=l) t 

c+\0) i = \S=l,M=-l) i 

4\0)i=\S=l,M = 0)i 
and |0)j is the particle vacuum on site i. 

We introduce the variational trial wavefunction 

|T^) = e F |0). (59) 

where 

p = h) E n t + h) E E + fa) E E E + m 

i\ ii 12 ii 12 is 

N is the total number of sites in the SO (4) representation, |0) the particle 
vacuum and 

Q t = a n a t + friK + lk c t + S H d t 

Developing the exponential in Eq. (|59|) leads to 

|T( Ar) ) = [l + F+iF 2 + iF 3 + ...]|0). (61) 

Since the creation operators {Qf} are fermions and since the system 
should be occupied by two particles by site i (1 particle on each site of a 
rung i) the sum of terms in Eq. (J6U|) and Eq. (J61)) is restricted to a finite 
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number of terms which all contain exactly N distinct operators {Q^ }. This 
sum is written in a fixed order so that %\ < % 2 < i% < ... < i^. The state 



lY^) can be written as 



N 



l T r) = El E ^(p;ni,n 2 ,...,n^)/ (1) / (2) .../ (iV )fi+..^|0). (62) 

p=l P' m,n 2 ,...njv 

where the sum contains only the terms which generate N operators {f^ 1- } 
and 

B N (p;n 1 ,n 2 ,...,n N ) = (-1) P 



n\\...n p \ 

with the constraint ri\ + 2n 2 + .... + Nn^ = N and P the number of permu- 
tations such that the operators {fif} appear in increasing order with respect 
to their index, i x < i 2 < 13 < ... < In- 

In a realistic description the correlation functions fm could depend on the 
relative distance between the particles on the different sites, = fn?'"' H \ 
We do not take this fact into account. It is sensible to believe that this 
point is not crucial for the aim which we fix here, i.e. to show the qual- 
itative behaviour of the system as already claimed before. It makes also 
physically sense to consider that the coefficients fm, (I > 2) which cor- 
relate / particles are a product of two-body terms f\i) = C(7,2)/( 2 ) where 
C(J,2)=J!/(2!(J-2)!). 

Under these conditions 



N ' N 

l T ™> = £^ E ^(p;n 1 ,...,n^) / -/;f-»n^(^)] n; l^ iV) )- 

p=l r m n\,n2,--.,n N 1=2 

(63) 

where n^..N) — n 2 + ... +un with the constraint n\ +2n 2 + ....Nun = N. 



8.2 Variational wavef unction. 

The structure of the correlated wavefunction |Tq^) is a product of a term 
which contains the correlations and the mean-field wavefunction |$q^). In its 
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whole generality the search of extrema of the energy 

is in principle given by varying the coefficients {Pi}, {7i}, {St} and 

/(!)) /(2) under the normalization constraint 



ri N) \r { N) ) = i. (64) 



One may consider a two-step procedure suggested by the structure of 
the wavefunction |Tq ) in Eq. (63). A minimization on the wavefunction 
1$^) is done first on E$ = (^ { N) \H^\^ { N) ) as before in section 4. Then 
Eq. (J64j) is verified if 

1 N 

C = ^^^n[C(/,2)]^(p;ni,... ) n iV )S^;n' 1) ...,n^)(65) 

P,p' P ' P ' rnX 1=2 

i i 
pni+rij , n (2..N)+n (2 N) _ 

J (1) J (2) ~ 1 

where the sums over {rii} are such that n\ + 2ri2 + .... + Nn^ = N and the 
same for {n^}. 

C should then be fixed by looking for values of /m and /( 2 ) such that the 
normalization constraint given by Eq. (j64*| is realized and is minimized. 

By construction it appears that all the reasoning concerning the nature 
of the spectrum which was made in sections 4 and 5 remains then valid as 
long as solutions exist. Hence the same conclusions as those drawn there can 
be taken up in the present case. 

If such a solution can be found it ensures the existence of a gap for any 
space dimension. 



9 Summary and conclusions. 

We introduced a variational approach in order to describe the low-lying en- 
ergy states of frustrated quantum spin systems. We started with spin-1/2 
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ladders and extended the approach to even-legged systems in any space di- 
mension d. 

In a first step the postulated trial wavef unction was written as a superpo- 
sition of spin singlet and triplet states in the framework of SO (4) symmetry. 
The variational ansatz limited the wavefunction to a mean-field product of 
site states formally similar to the ansatz used in high-T c theory jJT]. The 
local amplitudes of the site state components were correlated through near- 
est neighbour contributions. We looked for energy extrema and analyzed the 
conditions under which the variational procedure leads to energy minima. 
The realization of this property depends on the parameter subspace in which 
the coupling strengths of the spin-spin interactions are located. If the energy 
extrema do not correspond to minima the variational ansatz used for the trial 
wavefunction does not describe the structure of the physical ground state. In 
these cases the structure of this state must be different from the postulated 
ansatz. 

In the cases where the trial wavefunction describes a state of minimum 
energy the present analysis shows that spectra of spin- 1/2 quantum spin sys- 
tems with two legs are gapped. The approach can in principle be generalized 
to spin-spin interactions of any range and any type of frustration if the pos- 
tulated variational ground state wavefunction leads to an energy minimum. 
For specific values of the coupling constants the gap may vanish leading to 
degeneracies in energy. 

The extension of the formalism to even-legged systems of dimension d > 2 
was developed in section 5. The pertinence of such an extension which has 
been discussed in section 6 is warranted as long as the 5*0(4) symmetry 
which characterizes even-legged ladders remains a symmetry which governs 
the system at higher space dimensions. 

In a further step we suggested an extension of the variational trial wave- 
function to include explicit correlations beyond the mean-field approxima- 
tion. We introduced simplifications to the most general form of this wave- 
function which may have quantitative but not qualitative effects on the con- 
clusions. It is seen that the corresponding spectra may again show the char- 
acteristic gap obtained with the mean-field variational ansatz. 
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